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HIGHLIGHTS 


•  Proposed  a  unified  OCV  model  for  Li-ion  battery  SOC  estimation  and  SOH  monitoring. 

•  Considered  Li-ion  staging  phenomenon  during  intercalation  in  model  parametrization. 

•  The  model  improves  SOC  estimation  accuracy  compared  to  existing  OCV  models. 

•  The  model  captures  battery  electrochemical  properties  at  different  aging  stages. 

•  Established  relations  between  the  model  parameters  and  cell  performance  degradation. 
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Open-circuit-voltage  (OCV)  data  is  widely  used  for  characterizing  battery  properties  under  different 
conditions.  It  contains  important  information  that  can  help  to  identify  battery  state-of-charge  (SOC)  and 
state-of-health  (SOH).  While  various  OCV  models  have  been  developed  for  battery  SOC  estimation,  few 
have  been  designed  for  SOH  monitoring.  In  this  paper,  we  propose  a  unified  OCV  model  that  can  be 
applied  for  both  SOC  estimation  and  SOH  monitoring.  Improvements  in  SOC  estimation  using  the  new 
model  compared  to  other  existing  models  are  demonstrated.  Moreover,  it  is  shown  that  the  proposed 
OCV  model  can  be  used  to  perform  battery  SOH  monitoring  as  it  effectively  captures  aging  information 
based  on  incremental  capacity  analysis  (ICA).  Parametric  analysis  and  model  complexity  reduction  are 
also  addressed.  Experimental  data  is  used  to  illustrate  the  effectiveness  of  the  model  and  its  simplified 
version  in  the  application  context  of  SOC  estimation  and  SOH  monitoring. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  growing  concern  over  oil  shortage  and  environmental  is¬ 
sues  has  greatly  accelerated  the  development  of  alternative  power 
and  energy  solutions  to  displace  fossil  fuel  in  recent  years  [1—3],  As 
an  important  component  in  the  clean  energy  intuitive,  energy 
storage  system  is  a  critical  enabler  for  the  next  generation  power/ 
energy  technologies.  Among  the  various  choices  of  energy  storage 
technologies,  it  is  well  recognized  that  battery  systems  can  offer  a 
number  of  high-value  opportunities  given  their  advantages  in 
portability  and  energy  efficiency,  particularly  for  mobile 
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applications  [4],  Given  the  fact  that  batteries  have  been  used  to 
power  an  increasingly  diverse  range  of  applications,  from  cell 
phones  to  electric  vehicles,  more  reliable  and  efficient  battery 
management  systems  (BMSs)  need  to  be  developed  [5], 

BMS  has,  among  many  functions,  two  important  tasks,  namely 
the  state-of-charge  (SOC)  estimation  and  state-of-health  (SOH) 
determination  [6,7],  SOC  is  commonly  defined  as  “the  percentage 
of  the  maximum  possible  charge  that  is  present  inside  a 
rechargeable  battery",  and  the  estimation  of  SOC  serves  as  the  fuel 
gauge  for  batteries.  On  the  other  hand,  SOH  is  “a  ‘measure’  that 
reflects  the  general  condition  of  a  battery  and  its  ability  to  deliver 
the  specified  performance  in  comparison  with  a  fresh  battery"  [8], 
Typically,  the  quantitative  definition  of  SOH  is  based  either  on  the 
battery  capacity  or  the  internal  resistance,  depending  on  specific 
applications  [7],  Due  to  the  requirement  for  online  estimation  of 
both  SOC  and  SOH  during  system  operation,  the  development  of  an 
simple  but  accurate  battery  model  is  crucial  for  an  effective  and 
robust  BMS. 
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Many  studies  of  battery  SOC  estimation  have  been  reported  in 
the  literature  (see  Ref.  [9]  and  references  therein).  Most  of  those 
methods  are  model-based  approach  and  require  an  accurate  open- 
circuit-voltage  (OCV)  model  which  relates  OCV  to  SOC  [9-11].  The 
OCV— SOC  function  is  implemented  in  the  BMS  either  as  a  look-up 
table  or  an  analytical  expression,  while  the  latter  has  several  ad¬ 
vantages  including  computational  efficiency  (since  no  interpolation 
is  needed)  as  well  as  the  convenience  for  analysis.  Because  of  the 
wide  flat  region  on  the  OCV-SOC  curve  for  lithium  iron  phosphate 
batteries,  a  small  mismatch  in  OCV  fitting  may  cause  a  large  devi¬ 
ation  in  SOC  estimation.  Therefore,  improving  OCV  models  can 
contribute  significantly  to  increasing  the  accuracy  of  SOC 
estimation. 

In  Ref.  [12],  several  phenomenological  OCV  models,  which  are 
built  with  curve  fitting  without  considering  the  complex  battery 
physical  behavior  during  the  lithium-ion  intercalaction/dein- 
tercalation  process  [1,13],  are  summarized.  Although  those  models 
represent  part  of  the  nonlinear  characteristics  of  battery  OCV,  they 
fail  to  capture  the  voltage  plateaus  and  transitions  (please  see  Fig.  1 ) 
that  correspond  to  the  staging  phenomenon  at  the  graphite  anode 
[7,14-18], 

The  modeling  of  OCV-SOC  function  is  also  found  to  be  crucial 
when  it  comes  to  battery  SOH  monitoring,  as  OCV  data  often  reflect 
battery  aging  and  performance  degradation  [19],  In  particular,  by 
differentiating  the  battery  charged  capacity  (Q)  with  respect  to  the 
terminal  voltage  (V),  the  voltage  plateaus  on  the  OCV  curve  is 
transformed  into  clearly  identifiable  dQJdV  peaks  on  the  incre¬ 
mental  capacity  (IC)  curve  (see  Fig.  1)  [7,14,18],  This  so-called  in¬ 
cremental  capacity  analysis  (ICA)  technique  helps  to  detect  gradual 
changes  in  cell  behavior,  based  on  OCV  data  collected  from  life  cycle 
test,  with  greater  sensitivity  than  those  based  on  conventional 
methods  [5],  ICA  is  applicable  to  lithium  ion  batteries  that  use 
graphite  as  the  negative  electrode  material  since  it  originates  from 
the  study  of  the  lithium  intercalation  process  [14—17],  The  efficacy 
of  ICA  has  been  shown  with  various  lithium  ion  batteries  (LiFeP04, 
LiNMC,  LNCAO,  LiMn204,  etc.)  [7,18,20-23],  In  Ref.  [23],  the  aging 
signature  extracted  using  ICA  is  further  amplified  through  the  delta 
differential  capacity  analysis  technique,  which  enables  people  to 
detect  small  battery  degradation  in  a  shorter  time  interval 
[14,24,25],  One  alternative  way  of  extracting  battery  aging  signa¬ 
ture  and  analyzing  battery  SOFI  from  the  voltage  curve  is  the  dif¬ 
ferential  voltage  analysis  (DVA).  In  contrast  to  the  ICA,  the  DVA  is 


Measured  Voltage  (V) 


curve  at  1/20  C  rate. 


based  on  the  dV/d Q  curve.  The  accuracy  of  DVA  based  capacity 
fading  prediction  has  been  shown  in  Refs.  [26,27], 

One  major  difficulty  in  performing  ICA  and  DVA  is  the  sensitivity 
to  noise  in  battery  voltage  measurement  [7,21,23],  Since  all  the 
peaks  on  an  IC  curve  lie  within  the  flat  region  of  the  OCV  curve, 
computing  the  derivatives  directly  from  the  data  could  lead  to 
inaccurate  and  undesirable  results,  even  after  careful  data  filtering 
[7],  This  issue  could  be  addressed  by  implementing  specially 
designed  high  precision  testing  equipment  [23,25],  However,  for 
on-board  application,  the  measurement  precision  is  usually  limited 
and  appropriate  data  processing  is  required. 

By  constructing  an  OCV  model  that  represents  the  underlying 
physical  phenomenon  of  lithium-ion  intercalaction/dein- 
tercalation  process  and  therefore  captures  the  subtle  transitions 
over  the  flat  region,  one  can  take  advantage  of  the  analytical 
form  of  the  model  to  extract  useful  information  associated  with 
battery  electrochemical  properties  and  aging  status  using  ICA 
[5,20],  The  effect  of  measurement  noise  could  also  be  mitigated 
by  fitting  the  raw  data  with  a  specially  structured  parametric 
model.  Therefore,  for  both  SOC  estimation  and  SOH  monitoring 
through  ICA,  an  OCV  model  that  can  capture  the  physical  process 
over  the  flat  region  of  the  OCV-SOC  curve  and  be  capable  of 
supporting  ICA  is  needed. 

In  this  paper,  we  propose  a  unified  OCV-SOC  model  that  is 
intended  for  both  SOC  estimation  and  SOH  monitoring.  The  original 
idea  was  initially  proposed  in  Ref.  [37],  The  study  is  based  on  a 
battery  life  cycle  test  data  set  collected  from  eight  A123  APR18650 
cells,  which  use  LiFeP04  as  the  positive  electrode  material  and 
graphite  as  the  negative  electrode,  over  a  period  of  18  months 
[7,28],  The  test  data  used  for  this  study  are  acquired  through  a 
battery  test  bench,  which  includes  an  Arbin  BT2000  tester,  a  ther¬ 
mal  chamber  for  environment  control,  a  computer  for  user-ma¬ 
chine  interface  and  data  storage,  a  switch  board  for  cable 
connection,  and  battery  cells  [7,28],  The  data  acquisition  system  has 
a  logging  frequency  of  10  Hz,  and  the  measurement  precision  of 
both  current  and  voltage  is  0.02%  (i.e.,  1  mV  for  voltage 
measurement). 

The  remainder  of  this  paper  is  organized  as  follows.  We  first 
propose  the  development  and  parametrization  of  the  new  OCV 
model  in  Section  2.  Then  we  compare  the  new  OCV  model  with 
other  existing  models  in  terms  of  OCV  data  fitting  and  SOC  esti¬ 
mation  accuracy  in  Section  3.  In  Section  4,  we  elucidate  the  appli¬ 
cation  of  the  OCV  model  for  SOH  monitoring,  followed  by 
parametric  analysis  and  model  simplification  in  Section  S.  The 
conclusions  are  given  in  Section  7. 


Fig.  2.  Fitting  result 


Fig.  1.  ICA  performed  on  charging  voltage  i 


:  with  the  proposed  OCV  model. 


2.  A  new  open-circuit-voltage  model 

We  propose  a  new  OCV  model  structure  that  can  capture  the 
phenomenological  characteristics  associated  with  the  lithium-ion 
intercalaction/deintercalation  process  while  fitting  the  OCV— SOC 
data.  A  specific  model  structure  is  chosen  to  enforce  the  model  to 
exhibit  plateaus  on  the  OCV  curve.  These  plateaus  are  results  from 
the  lithium-ion  staging  phenomena,  and  have  been  irrefutably 
observed  from  experimental  data.  The  parameters  of  the  proposed 
model  will  be  identified  with  experiment  data.  Instead  of 
measuring  the  OCV  after  a  long  relaxation  period  at  different  SOC 
levels  that  span  the  entire  range,  the  data  used  are  collected  by 
charging/discharging  a  LiFeP04  battery  cell  at  low  rate  of  1  / 20  C 
under  room  temperature.  The  voltage  data  obtained  by  this  slow 
charging/discharging  process  reflects  the  OCV  at  a  close-to- 
equilibrium  status  [18],  Although  the  close-to-equilibrium  OCV 
curve  shown  in  Fig.  2  may  be  affected  by  hysteresis  and  diffusion 
voltage  due  to  lack  of  relaxation,  it  can  sufficiently  represent  the 
generic  electrochemical  properties  for  ICA  and  quantifying  capacity 
fading  as  presented  in  Ref.  [  18  ].  The  upper  and  lower  cut-off  voltage 
limits  are  set  to  be  3.6  V  and  2.1  V  respectively.  The  total  charge 
stored  between  the  upper  and  lower  limits  is  defined  as  the  battery 
capacity. 


OCV(z)=/f0  +  /<iTT-^ 


z\+ea2(z-P2) 

1  +  I(5Z 


(1) 


where  z  is  the  SOC,  I(0~  5  are  the  linear  parameters,  and  ai~4  and 
~2  are  the  nonlinear  parameters.  This  parametrization  shares 
features  with  some  of  the  electrochemical  model  of  LiFeP04  pro¬ 
posed  in  Ref.  [1  ].  More  specifically,  our  model  can  be  interpreted  as 
a  simplified  composition  and  generalization  of  the  cathode  and 
anode  OCV  functions  in  Ref.  [1]  (where  hyperbolic  and  exponential 
functions  are  used). 

An  example  of  a  complete  data  set  for  model  identification  is 
shown  in  Fig.  2,  where  the  collected  charging/discharging  data  are 
averaged  to  yield  OCV  data.  The  parameters  in  Eq.  (1)  are  deter¬ 
mined  using  the  Matlab  curve  fitting  toolbox.  We  focus  on  fitting 
the  data  in  the  middle  SOC  range,  namely  10%-90%  SOC,  given  that 
in  practice  only  data  in  the  middle  SOC  range  are  available.  More¬ 
over,  this  allows  us  to  achieve  higher  accuracy  in  the  range  where 
the  model  will  be  most  useful.  The  fitted  OCV— SOC  curve  is  plotted 
in  Fig.  3,  with  a  reference  made  to  the  averaged  measured  OCV  data. 
The  rms  and  maximum  values  of  fitting  error  is  1  mV  and  2.5  mV, 
reduced  from  2.3  mV  (rms)  and  7  mV  (maximum)  that  were  ob¬ 
tained  when  the  full  SOC  range  is  used.  Table.  1 


2.1.  Open-circuit-voltage  parametric  model 

Based  on  the  studies  presented  in  Refs.  [7,14,18],  there  are 
generally  three  observable  voltage  plateaus  and  two  transitions, 
which  are  important  in  characterizing  the  OCV  and  battery  capac¬ 
ity,  over  the  flat  area  on  the  OCV  curve  of  LiFeP04  battery  cells.  In 
our  model,  those  plateaus  and  transitions  are  represented  by  the 
sigmoid  functions  as  follows, 


2.2.  Model  comparison  and  analysis 

The  new  OCV  model  is  compared  in  this  section  with  the  five 
different  models  summarized  in  Ref.  [12],  where  the  polynomial 
model  is  reported  to  be  the  most  accurate.  All  the  parameters  in 
those  five  OCV  models  are  refitted  for  the  data  presented  in  Fig.  2 
using  Matlab  curve  fitting  toolbox,  the  results  with  their  rms  and 
maximum  errors  are  shown  and  summarized  in  Fig.  4  and  Table  2.  It 


Table  1 

OCV  models  summarized  in  Ref.  [12]  together  with  the  new  model. 


# 

OCV  models 

Ref. 

1 

OCV(z)  =  K0-  K\lz  -K2Z  +  if3ln(z)  +  JCJn(l  -  z) 

[32] 

2 

OCV(z)  =  Ko+KjO  -  e~“'z)  +  K2  (l  -  e-S)  +  K3z 

[33] 

3 

OCV(z)  =  K0  +  K^e-“'d~z)  -K2/z 

[34] 

4 

OCV(z)  =  K0  +  Ki  e-“'z  +  K2z  +  K3z2+  JQz3 

[35] 

5 

OCV(z)  =  K0  +  KjZ  +  K2Z2  +  K3z3  +  K4Z4  +  K5Z5  +  KeZ6 

[36] 

6 

OCV(z)  =  Kq  +Kjj+eaJwy  +*k1+e4n2)  +  ki  i+e°U  ')  +*4i+ev  + 

Proposed  model 
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Max  i 


12.6 

8.9 

21.3 

12.7 

73 

2.5 


(mV) 


OCV(V) 


models  #2  and  #5,  dQJdV  is  not  defined  at  3.35  V  because  of  a 
singularity  in  the  math  expression.  Consequently,  spikes  in  both 
directions  are  shown  in  Fig.  5  for  the  two  IC  curves,  corresponding 
to  models  #2  and  #5. 


3.  State-of-charge  estimation  based  on  extended  Kalman 
filter 


The  extended  Kalman  filter  (EKF)  based  approach  discussed  in 
Ref.  [29]  is  used  to  illustrate  the  implementation  of  the  new  OCV 
model  for  SOC  estimation.  The  battery  test  data  for  this  study  is 
collected  through  the  experimental  set-up  introduced  in  Ref.  [28] 
and  the  first-order  RC  model  [30]  is  used  for  calculating  the  SOC 
with  the  following  algorithm: 


RC  Model : 

zk+l  -  f(Zk,k)  +  Wk , 

Vk  =  g(zk,Ik)  +  vk, 

EKF  Algorithm : 

a  9fcw*-»y| 

*fc-l  ~  WXm* 

Cfc  =  ?g^M|  ^  , 

k  'Zk-l=Zk\k-l 

zk|k-i  =f(zk-  i,4-i). 

Pfc|k-i  =  Ak  ]Pk  iA[  ,  +  Q. 


(2) 


Fig.  5.  Comparison  of  IC  curves  from  fitted  OCV. 


should  be  noted  that  none  of  the  five  OCV  models  is  suitable  for  IC 
analysis,  as  their  model  structures  do  not  take  into  account  the 
staging  mechanism  in  intercalation/deintercalation  process.  One 
can  see  that  the  new  OCV  model  proposed  in  Eq.  (1)  has  better 
fitting  accuracy  than  all  those  five  models.  Consequently, 
improvement  in  the  SOC  estimation  results  may  also  be  expected 
when  the  new  OCV  model  is  incorporated. 

The  IC  curves  based  on  different  OCV  models  are  shown  in  Fig.  5. 
It  can  be  observed  that  the  new  model  (model  #6)  captures  the  two 
IC  peaks  associated  with  staging,  whereas  all  other  models  only 
show  one  peak  on  the  IC  curves.  Therefore,  model  #6  can  extract 
more  capacity  information  from  the  OCV  data  compared  to  other 
models  and  will  be  used  for  further  analysis  in  SOH  monitoring.  For 


0  2000  4000  6000  8000  10000 


Gk  =  Pk\k-iCTk[ckPk]k_,CTk+v]  \ 
zk  =  zk|k-i  +  Gk  [Vk  -  g(zk |k_i ,  4)  ] , 

Pfc  =  [/  -  GkCk]Pklk  , , 
where 

/(zfc>  4!  *  zk  — 

g(zk>  4)  =  OCV(Zfc)  +  Vrck  +  R0Ik, 

Kc,k+ 1  =  exp (r,Ac)  Vrck  +  fit  [1  -  exp(Rfct)]4- 

zk  is  the  SOC,  rji  is  the  charging/discharing  efficiency,  At  is  the 
time  step,  Qc  is  the  battery  capacity,  Ik  is  the  input  current,  Vk  is  the 
total  terminal  voltage,  Vrc_k  is  voltage  of  the  RC  circuit,  R0.  Ri  and  C 
are  battery  internal  resistance  and  capacitance,  Q  is  the  error 
covariance  of  process  noises,  and  V  is  the  error  covariance  of 
observation  noises  [28],  The  OCV  function  is  represented  by  the 
model  given  by  Eq.  (1). 

In  this  study,  all  the  parameters  in  the  RC  model  are  set  to  be 
constants  without  performing  iterative  optimizations  as  discussed 
in  Ref.  [12], 

The  battery  data  tested  under  the  Federal  Urban  Driving 
Schedule  (FUDS)  is  shown  in  Fig.  6.  As  in  Ref.  [29],  the  SOC  data 
obtained  by  Coulomb  counting  in  the  Arbin  battery  test  bench, 
where  high  precision  current  sensor  is  used,  is  assumed  to  be  the 
“true"  SOC  and  used  as  a  reference  for  performance  evaluation. 
Therefore  the  SOC  estimation  error  is  defined  as, 

esoc  ,k  =  zk“zk  (3) 

where  esoc,k  is  the  SOC  estimation  error,  zk  is  the  SOC  estimated  by 
the  EKF  algorithm  and  zk  is  the  reference  SOC  calculated  from  the 
measured  data.  Please  note  that  the  SOC  from  Arbin  Coulomb 
counting,  while  representing  the  best  reference  that  we  can  get,  is 


Fig.  6.  Battery  data  tested  under  FUDS. 
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Table  3 

Parameters  of  models#5  and#6  used  in  SOC  estimation. 
Model  parameters  Model#5 


Ko 

Ki 

K2 

k3 

k4 

k5 


3.0896 

1.1627 

-2.3821 

2.1870 

-0.5444 

-0.1939 

0.0582 

N/A 

N/A 

N/A 

N/A 


Model#6 


3.4002 

0.0080 

0.0785 

-0.2150 

-1.3032 

0.0891 

N/A 


40 


Parameters  of  EKF  based  SOC  estimator. 

Parameters 

Values 

Error  covariance  of  process  noises,  Q 

Error  covariance  of  observation  noises,  V 

(°o15  ?) 
0.01 

only  approximately  accurate  due  to  the  unavoidable  error  and  in¬ 
tegral  accumulation  of  the  error  in  the  current  measurement  [29]. 

The  model  parameters  used  in  the  EKF  based  SOC  estimation  are 
shown  in  Table  3. 

A  detailed  description  of  the  EKF  algorithm  can  be  found  in  both 
Refs.  [29]  and  [12],  The  parameters  of  the  EKF  (i.e.,  V  and  Qused  in 
Eq.  (2))  are  calibrated  using  the  data  collected  in  our  aging  tests, 
and  the  same  parameters  (listed  in  Table  4)  are  used  in  the  per¬ 
formance  evaluation  for  the  new  model  (#6)  and  the  polynomial 
model  (#5). 

Fig.  7  displays  the  EKF  based  SOC  estimation  results.  The  two 
plots  represent  two  cases  with  different  initial  SOC  error  (+10% 
and  -10%  respectively).  One  can  see  that  the  estimated  SOC 
converge  into  the  5%  estimation  error  bound  when  either  model  #5 
or  #6  is  used.  In  particular,  with  positive  initial  error,  the  SOC 
estimated  with  model  #6  converges  much  faster  than  the  SOC 
estimated  with  model  #5.  This  difference  in  convergence  rate  with 
positive  initial  error  can  be  explained  by  the  curve  fitting  results 
shown  in  Figs.  3  and  4,  where  model  #6  have  better  accuracy  than 
#5  in  the  high  SOC  region.  On  the  other  hand,  both  models  have 
approximately  the  same  accuracy  when  SOC  is  below  70%  and 
therefore  have  similar  convergence  rate  with  negative  initial  error. 


Fig.  7.  SOC  estimation  with  OCV  model  #6. 
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Fig.  8.  V-SOC  and  IC  curve  from  0.5C  charging  data. 


The  5%  error  bound  can  be  achieved  with  other  initial  error  (e.g., 
20%  or  30%). 

In  summary,  the  new  OCV  model  performs  well  in  EKF  based 
SOC  estimation.  The  model  delivers  accurate  estimation  results 
with  less  than  5%  errors  without  extensive  calibration  and  training. 
The  results  presented  in  this  section  demonstrates  the  potential  of 
implementing  this  new  OCV  model  for  SOC  estimation  problems. 

4.  State-of-health  monitoring  based  on  incremental  capacity 
analysis 

Even  though  ICA  was  originally  proposed  for  “close-to-equilib- 
rium”  conditions,  it  was  shown  in  Ref.  [7]  that,  for  normal  charging 
data,  the  peaks  on  the  IC  curve  can  also  be  identified  and  they 
reveal  significant  information  about  battery  SOH.  Hence,  the  results 
presented  in  this  section  for  ICA  are  all  based  on  0.5C  battery 
charging  data  (same  as  in  Ref.  [7]).  It  will  be  shown  that  the  pro¬ 
posed  parametrization  is  applicable  to  both  OCV  data  and  normal 
charging  data  (i.e.,  charged  at  0.5C). 

As  shown  in  Fig.  8,  where  the  charging  data  is  fitted  with  the 
model  given  by  Eq.  ( 1 )  with  OCV  replaced  by  the  terminal  voltage 
collected  at  normal  charging  rate,  we  can  obtain  the  IC  curve 
directly  by  differentiating  the  analytical  V-SOC  function  and  then 
taking  the  reciprocal. 

4.1.  Incremental  capacity  analysis  at  different  aging  stages 

The  IC  curves  of  battery  charging  data  collected  at  different 
aging  cycles  under  35  °C  environment  condition  are  shown  in  Fig.  9. 
The  numerical  values  in  the  plot  legends  represent  the  aging  cycle 
number.  The  1C  peaks  highlighted  by  the  circle  in  Fig.  9  clearly 
shows  a  monotonic  decreasing  trend  as  battery  ages,  which  implies 
a  battery  degradation  that  is  mainly  related  to  the  loss  of  active 
material  at  the  graphite  anode  [14,18], 

The  IC  curves  from  four  different  cells  are  plotted  in  Fig.  10, 
where  the  plots  focus  on  the  second  peak  on  the  IC  curve  and  depict 
a  consistent  trend  in  capacity  fading  as  reflected  by  the  decreasing 
IC  peaks. 

To  obtain  the  correlation  between  the  cell  capacity  fading  and  1C 
peak,  the  peak  values  of  eight  cells  at  different  aging  stages  are 
identified.  The  normalized  peak  values  of  one  of  the  cells  (red  dots 
in  the  web  version  in  Fig.  11)  is  correlated  with  the  faded  cell  ca¬ 
pacities,  and  used  to  construct  the  battery  capacity  estimation 
curve  (red  dashed  curve  in  Fig.  11 )  as  done  in  Ref.  [7],  The  data  from 
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Fig.  9.  1C  curves  of  data  at  different  aging  stages  at  35  “C.  Fig. 11.  Correlation  between  faded  battery  capacity  and  IC  peak  values  at  35  °C 


other  cells  are  also  plotted  to  verify  the  prediction  capability  of  the 
estimation  curve.  The  result  shows  the  new  parametric  model  in 
the  form  of  Eq.  ( 1 )  can  be  used  for  capacity  based  battery  SOH 
monitoring  with  good  accuracy. 

4.2.  Incremental  capacity  analysis  at  different  temperature 


second  peak  shows  the  aging  status.  It  suggests  that  the  test  data 
corresponding  to  higher  temperature  will  be  more  reliable  for  ICA 
and  capacity  fading  identification.  Therefore,  even  though  the 
temperature  dependence  is  not  parametrized  explicitly,  the  pro¬ 
posed  model  is  capable  of  capturing  the  influence  of  temperature 
difference  by  adapting  all  parameters. 


Test  data  were  also  collected  at  different  temperature 
throughout  the  battery  aging  process.  ICA  was  performed  with 
respect  to  different  temperatures  to  evaluate  and  confirm  the 
sensitivity  to  temperature  and  results  are  presented  in  Fig.  12.  The 
intensities  of  the  IC  peaks  reduce  as  temperature  drops,  which 
could  be  attributed  to  the  slow  lithium-ion  diffusion  or  lithium 
plating  phenomenon  at  low  temperature  as  elucidated  in  Ref.  [31], 
When  the  temperature  is  too  low,  say  at  10  °C,  the  second  peak 
could  disappear  from  the  IC  curve.  However,  since  the  operating 
temperature  of  battery  systems  are  usually  maintained  through 
integrated  thermal  control  in  the  BMS,  the  disappearance  of  IC 
peaks  at  low  temperature  should  not  be  a  concern. 

Fig.  13  presents  the  IC  curves  at  different  temperature  at 
different  aging  stages.  A  temperature  dependence  of  the  IC  peaks 
and  the  model  parameters  may  be  evaluated.  The  first  peak  on 
these  four  plots  shows  consistent  temperature  sensitivity  while  the 
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cells  at  35  °C. 


5.  Other  signatures  on  the  IC  curve  and  their  correlation  with 
capacity  fading 

Other  than  using  the  IC  peak  value  decreasing  as  the  indication 
of  battery  capacity  fading,  we  also  explore  a  different  approach  that 
uses  the  IC  peak  locations  to  evaluate  the  degradation.  Since  the  IC 
peaks  represent  the  phase  transitions  during  the  intercalaction/ 
deintercalation  process  and  this  electrochemical  reaction  is 
affected  by  battery  capacity  loss,  the  locations  of  the  peaks  might 
also  have  a  correlation  with  the  faded  capacity  and  can  give  a  good 
estimation  of  the  loss. 

Consider  the  IC  as  a  function  of  SOC  (rather  than  measured 
voltage),  as  shown  in  Fig.  14.  Note  that  the  first  IC  peak  stays  at 
approximately  the  same  location  whereas  the  second  peak  is 
moving  towards  the  left  as  the  battery  ages  and  capacity  decreases. 
By  normalizing  both  the  second  IC  peak  location  and  battery  ca¬ 
pacity,  we  obtain  the  correlation  between  the  two  values  in  Fig.  15 
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Fig.  12.  IC  curves 
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Fig.  13.  IC  curves  of  data  at  different  temperature  and  ages. 


(data  from  different  cells  are  color  coded  differently  (in  the  web 
version)),  which  confirms  that  the  location  of  the  second  IC  peak 
could  also  be  a  indication  of  battery  capacity  fading. 

However,  in  order  to  accurately  measure  this  peak  location  in 
terms  of  charged  capacity,  it  is  required  to  start  the  battery  charging 
from  zero  SOC,  which  is  not  always  practical  in  real  operation.  This 
problem  can  be  addressed  by  using  the  SOC  span  between  the  two 
consecutive  1C  peaks, 


The  analysis  suggests  that  the  signatures  on  the  IC  curves, 
namely  the  second  IC  peak  and  the  distance  between  two  peaks, 
can  be  used  as  indicators  for  battery  aging.  The  two  correlations 
based  on  IC  peak  locations  are  not  as  strong  and  robust  compared  to 
the  results  based  on  peak  values  shown  in  Fig.  11.  Nevertheless,  we 
believe  that  those  correlations  can  provide  useful  information  for 
battery  SOH  monitoring  from  different  perspectives.  It  is  also 
possible  to  improve  the  peak  location  based  results  when  more 
accurate  current  sensors  are  available. 


d,c 


Pi -Pi 


(4) 


6.  Parametric  analysis  and  model  simplification 


where  die  is  the  SOC  interval  between  two  consecutive  peaks,  Pi 
and  P2  are  the  locations  of  the  two  IC  peaks  respectively.  The  value 
of  die  is  more  easily  obtainable  as  long  as  the  charging  starts  below 
the  SOC  value  that  corresponds  to  the  first  IC  peak.  Typically,  the 
first  peak  on  the  IC  curve  appears  around  SOC  =  0.45.  Therefore,  for 
all  the  data  we  have,  this  condition  is  satisfied  as  long  as  the 
charging  covers  40%— 85%  SOC.  The  plot  with  normalized  dIC  and 
battery  capacity  is  shown  in  Fig.  16.  This  correlation  obtained  using 
die  is  looser  than  using  P2  only,  as  the  values  of  both  P]  and  P2  are 
affected  by  measurement  noises  and  the  value  of  die  is  more  sus¬ 
ceptible  to  noise  and  modeling  error. 


Although  the  identification  results  using  the  new  OCV  model 
has  shown  its  effectiveness,  the  high  number  of  parameters  and 
nonlinearity  of  the  model  may  present  a  challenge  for  on-board 
implementation.  A  parametric  analysis  is  performed  to  charac¬ 
terize  the  correlations  between  model  parameters  and  battery  ag¬ 
ing,  evaluate  the  sensitivity  and  robustness  of  the  model  structure, 
and  better  understand  the  utility  of  the  proposed  model.  Moreover, 
the  analysis  results  may  help  to  identify  the  significance  of  each 
parameter  and  simplify  the  model,  so  that  the  computational  effi¬ 
ciency  in  estimation  can  be  improved. 

Consider  the  following  function,  which  is  derived  by  differen¬ 
tiating  OCV  with  respect  to  the  battery  SOC, 


i 

! 
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Fig.  14.  IC  curve  plotted  versus  charged  capacity  data  at  35  "C 


Fig.  16.  Correlation  between  IC  peak  location  difference  and  battery  faded  capacities. 
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G(z)  =  ^OCV(z) 


(1  +  e“’(z-W) 


(1  +  e“3(z-i))‘ 


(1  +  ea2(z-&))2 

g4e“4(^fa) 


Note  that  the  IC  curve  is  the  inverse  of  G(z),  the  locations  and 
values  of  IC  peaks  may  be  calculated  by  letting  the  derivative  of  G(z) 
equal  to  zero  and  solving  the  resulting  algebraic  equation.  The 
resulting  equation,  however,  is  highly  nonlinear  and  cannot  be 
easily  solved  analytically  in  general.  Here,  an  alternative  approach 
that  separates  and  analyzes  the  terms  in  the  Eq.  (5)  individually  is 
used  to  correlate  the  model  parameters  and  the  solution  of  G(z). 

To  understand  the  contribution  of  each  term  of  (5),  we  first 
define  the  following, 


/l(z)  = 

(l+e“uz 

-*)) 

h(z)  = 

K  ^ 

-h) 

(nW* 

h(z)  = 

-vs 

fM  ^ 

-K‘iSS 

Through  analytical  and  numerical  analysis, 
following  observations: 


(6) 


have  the 


/2  and/4  decay  to  zero  at  the  high  SOC  range,  where  the  IC  peak  of 
interest  occurs.  Therefore,  they  do  not  contribute  to  the  IC  peaks 
and  their  locations.  Our  analysis  will  be  focused  on/i  and/3: 

The  parameter  K5  is  a  constant  and  therefore  it  does  not  need  to 
be  considered  here; 

Numerical  results  show  that  I( \  and  K3  do  not  vary  much 
compared  to  a\  and  a3\ 

Our  analysis  show  that  the  IC  peak  has  about  equal  sensitivity  to 
a\  and  a3.  The  sensitivity  analysis  is  performed  by  numerically 
evaluating  the  partial  derivative  of  the  IC  function,  defined  by  1  / 
C(z)  with  G(z)  given  by  (5),  with  respect  to  the  parameters  a\ 
and  a3  at  the  z  value  where  the  second  IC  peak  is  located. 


With  the  above  observations,  we  identified  a  simplified  model 
in  the  following  form: 


■  e*  (*-*)' 


1+e-Hy 


where  the  model  structure  contains  only  4  changing  parameters 
(i.e.,  Ki,  /</  a\  and  0:3),  and  other  parameters  are  kept  constants  as 
batteries  age. 

In  Fig.  17,  one  can  see  that  the  IC  curves  obtained  using  the 
simplified  OCV  model  show  qualitatively  the  same  trend  as  the 
ones  obtained  by  the  full  model  (see  Fig.  9)  when  the  battery  ages. 
Table  5  shows  the  comparison  of  fitting  accuracy  between  the  full 
and  the  simplified  model.  In  the  worst  case,  the  rms  error  is  less 
than  2  mV.  Therefore,  the  simplification  preserves  the  high  fitting 
accuracy  provided  by  the  full  model. 

The  parameters  for  the  simplified  model  are  listed  in  Table  6, 
where  one  can  see  the  variations  in  ai  and  a3  agree  with  our 
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Fig.  17.  IC  curves  obtained  using  simplified  OCV  model  at  35  °C. 


Table  5 

Comparison  of  fitting  accuracy  between  full  and  simplified  model  using  RMS  error. 
Cycle#  Full  model  (mV)  Simplified  model  (mV) 

1000  0.760  0.805 

1400  0.925  0.980 

1900  0.946  1.247 

2300  1.024  1.559 


parametric  analysis.  Moreover,  it  suggests  that  the  values  of  the 
deteriorating  peaks  in  the  IC  curves  might  be  estimated  directly 
from  the  parameters  with  the  following  relationship, 

q  =  8q  +  6\Oi\  +  $2^1  +  63013,  (8) 

where  q,  defined  as  q  =  p  -  polpo,  represents  the  amount  of 
degradation,  p  represents  the  IC  peak  values  and  po  is  the  initial 
peak  value,  do- 3  ate  the  identification  parameters  in  this  equation. 
K3  is  not  included  in  this  equation  as  it  is  insensitive  to  the  variation 
of  IC  peak  values.  Same  as  our  analysis  in  the  Subsection  4.2,  the 
values  of  parameters  do -3  depend  on  the  operating  temperature. 
For  on-board  battery  SOH  diagnosis,  thermal  control  should  be 
integrated  as  part  of  the  BMS.  Therefore,  consistent  temperature 
could  be  expected  when  the  SOH  monitoring  is  performed. 

Given  a  set  of  n  data  points,  pu  an,  Ku,  a3 i  =  1  ,...,n,  from  one 
battery  cell  at  different  aging  stages,  the  values  of  do~3  can  be  easily 
calculated  from  the  least  squares  method, 

6  =  (V<J>)~Vq, 

where 

e  =  [60,63,82,63],  (9) 

®  =  2,-,0„]T, 

q  =  §fi,g2,  •••»«#!* 

h  =  il,«H.Kn,«3t]T. 


Table  6 

Parameters  of  the  simplified  OCV  model. 


Cycle#  a,  K,  a3  K3 

1000  -12.57  0.0446  25.9352  -0.2459 

1400  -12.29  0.0467  23.2749  -0.2435 

1900  -11.57  0.0497  20.4778  -0.2413 

2300  -11.4  0.0531  18.5991  -0.2447 
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Fig.  20.  Correlation  between  capacity  fading  and  IC  peak  reduction. 


The  identification  results  of  (8)  are  very  accurate  as  shown  in 
Fig.  18,  where  the  battery  testing  data  of  one  single  cell  (denoted  as 
cell  A)  at  all  different  ages  are  included. 

To  use  the  model  (8)  identified  from  cell  A  to  7  other  cells,  one 
can  use  ~3  directly  while  adjusting  do  such  that  the  value  of  q 
given  by  (8)  at  the  initial  stage  is  zero.  Fig.  19  show  that  the  data 
from  two  other  cells,  as  an  example,  can  also  be  fitted  accurately  by 
applying  d\  ~3  without  further  identification  effort.  In  other  words, 
the  IC  peak  values  of  different  cells  can  be  estimated/predicted 
based  on  the  parameters  $i~3  identified  from  one  reference  cell. 
Therefore  it  is  possible  to  perform  battery  SOH  monitoring  directly 
using  (8)  without  numerically  deriving  the  IC  curves. 

As  shown  in  Fig.  20,  the  results  from  Figs.  11  and  18  can  be 
combined  so  that  the  battery  capacity  fading  can  be  directly  pre¬ 
dicted  using  the  parameters  of  the  simplified  OCV  model.  One  can 
see  that  50%  reduction  in  1C  peak  value  to  equivalent  to  10%  battery 
capacity  loss. 

To  summarize,  the  simplified  OCV  model  is  capable  of  capturing 
the  battery  aging  features  with  fewer  parameters,  thereby 
improving  the  identification  efficiency  and  making  the  model 
implementable  in  on-board  BMS.  With  the  established  linear 
relationship  between  the  IC  peak  values  and  the  OCV  model  pa¬ 
rameters  in  (8),  the  battery  aging  information  can  be  revealed  by 
only  looking  at  the  parameters. 


Fig.  19.  Verification  of  (8)  using  parameters  0j  _3  identified  from  Fig.  18  with  data  from 


7.  Conclusions 

In  this  paper,  we  propose  a  new  OCV  parametric  model.  The 
special  parametrization  considers  the  staging  phenomenon  during 
the  lithium  intercalation/deintercalation  process  for  lithium-ion 
batteries,  therefore  leading  to  a  much  better  fitting  accuracy 
when  applied  to  experimental  battery  OCV  data.  The  new  para¬ 
metric  model  is  applied,  together  with  an  EKF,  to  SOC  estimation 
and  its  effectiveness  is  demonstrated  on  FUDS  cycle  data.  We  also 
show  that  the  parametric  model  can  be  applied  for  ICA  based  SOFI 
monitoring.  Through  comparisons  of  IC  curves,  the  model  is  shown 
to  be  capable  of  reflecting  battery  cell’s  electrochemical  properties 
at  different  operating  temperature  conditions  and  aging  stages. 
Moreover,  through  parametric  analysis,  we  are  able  to  simplify  the 
model  structure  and  develop  a  function  that  characterizes  the  re¬ 
lations  between  the  model  parameters  and  the  degradation  of 
battery  performance,  which  enables  us  to  monitor  battery  SOH  and 
estimate  capacity  fading  only  based  on  the  model  parameters. 
Future  research  will  focus  on  developing  more  efficient  algorithms 
for  the  determination  of  the  OCV  model  parameters. 
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